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Abstract 

The study of the two-shell system started in “Pair of null gravitating 
shells I and II” is continued. The pull back of the Liouville form to the con¬ 
straint surface, which contains complete information about the Poisson brack¬ 
ets of Dirac observables, is computed in the singular double-null Eddington- 
Finkelstein (DNEF) gauge. The resulting formula shows that the variables 
conjugate to the Schwarzschild masses of the intershell spacetimes are simple 
combinations of the values of the DNEF coordinates on these spacetimes at 
the shells. The formula is valid for any number of in- and out-going shells. 
After applying it to the two-shell system, the symplectic form is calculated for 
each component of the physical phase space; regular coordinates are found, 
defining it as a symplectic manifold. The symplectic transformation between 
the initial and final values of observables for the shell-crossing case is written 
down. 


1 Introduction 


The present paper is the third in a series devoted to the two-shell system. The 
first paper, Ref. [Q], and the second, Ref. P, will be referred to as I and II. In I, 
all classical solntions that have a regnlar center on the left are described, and the 
space of solutions is parametrized by three discrete and four continuous parame¬ 
ters. The space of solutions is a candidate for the physical phase space and the 
parameters are candidates for Dirac observables. In II, the action functional for a 
single shell due to Louko, Whiting and Friedman is generalized to any number 
of in- and out-going shells. The pull back ©p of the corresponding Liouville form 
0 to the constraint surface T is transformed into coordinates consisting of embed¬ 
dings, embedding momenta, and Dirac observables. Some general properties of the 
pull back, such as gauge invariance, are shown. They enable us to accomplish the 
transformation explicitly, and this is the main task of the present paper. 

The central idea of this paper is to calculate ©r in double-null Eddington- 
Finkelstein (DNEF) coordinates. This is a kind of gauge, but a singular one: Both 
the metric and the embeddings corresponding to this gauge are discontinuous at the 
shells and diverging at any Schwarzschild horizon. We shall show that this singu¬ 
larity does not influence the calculated value of ©p. That follows from the way ©p 
transforms under ordinary gauge transformations. What is the motivation for using 
the singular gauge? 

In fact, we first calculated it in regular gauges. The results revealed that the 
DNEF coordinates had a special role: their values at the shells appeared in the hnal 
formula for ©p. Moreover, the calculations were very long while the final formula 
was very simple, indicating that working with DNEF coordinates from the start 
could lead to simplihcations. This turns out to be right, although the calculation 
is still far from being trivial. One can perhaps say that the Liouville form chooses 
itself the spacetime coordinates in which it likes to be expressed. 

Originally, ©p has a form of an integral over a Cauchy surface. Each Cauchy 
surface has a boundary consisting of the regular center and the inhnity; each such 
surface intersects the shells at intersection points. Due to the discontinuity of the 
DNEF gauge at the shells, the contribution of each intersection point to 0p is non¬ 
zero, whereas the contributions from the regular center and the infinity in this gauge 
both vanish (cf. 0, where the only non-zero contribution comes from the infinity). 
Thus, ©p can be cast as a sum over all intersection points in which the summands 
have a standard form. In this way, a general formula can be shown to be valid for 
any number of in- or out-going shells. 

The sum over the intersection points is of course associated with a particular 
Cauchy surface. Let us consider two such surfaces. If the shells do not cross between 
these two Cauchy surfaces, then the corresponding summands have the same value on 
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each surface. However, if the shells do cross between the surfaces, the corresponding 
summands for the Cauchy surface below and above the crossing point are related 
by a highly non-trivial canonical transformation. 

The plan of the paper is as follows. In Sec. ^ the calculation of ©r in the DNEF 
gauge is justihed and accomplished. The result is a general formula that expresses 
©r in terms of some Dirac observables. These are the Schwarzschild masses of the 
intershell spacetime pieces and some combinations of the values of DNEF coordinates 
on these pieces at the shells. The canonical transformation between the observables 
below and above the crossing point is calculated. 

In Sec. I], the formula is specialized to the system of two shells. Some complete 
sets of Dirac observables as coordinates on the various components and regions of 
the physical phase space are considered. The physical phase space is then given the 
structure of a symplectic manifold. On each component of the physical phase space, 
we hnd a global chart with respect to which the components of the symplectic form 
are C°° and regular, in the sense that the matrix of the components has a nowhere 
vanishing determinant. In particular, in the case of shell-crossing, the “singular” 
case in which both shells lie on horizons (denoted by Coo in I) turns out to be a 
smooth surface in the phase space. Finally, Sec. ^ contains some conclusions and 
outlook, speculating about the prospective quantum theory. 

2 Calculation of the Liouville form in a singular 
gauge 

Let E be a Cauchy surface dehned by an embedding {U{p),V{p)). In II, the pull 
back ©r to the constraint surface T of the Liouville form © has been written as a sum 
of various contributions from different parts of E: First, there are contributions from 
each spacetime point (denoted by p) where a shell intersects E. If the intersection is 
with a single (in-going or out-going) shell, each point p contributes by a single term 
(cf. II, Eq. (46)) 

pdr . 

If the intersection is a point p where an in-going and an out-going shells cross each 
other, then the contribution is 


Pout dr -F pindr . 


Second, there are contributions from each connected volume cut out from E by the 
shell intersections. Each such volume contributes by the boundary terms 

{fdU -|- gdV + hidd' — ipdb)p=b — (fdU -|- gdV -|- hido^ — ipda)p=a , 
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where p = a and p = h, a < h, are the boundary points of the volume (cf. Eq. (60) 
of II), o* are Dirac observables, /, g, hi and ip are functions of o\ U{p) and V(p) 
defined by Eqs. (51)-(55) of II; 

/ = +mCo‘) . (1) 

9 = (2) 

hi = ^\n(-'h^+Hi(U,V,o‘) , (3) 

•p = + (4) 

-FU'-GV+ i>{U,V,o‘) , (5) 

The boundaries a and b can correspond either to the regular center p = 0 of S, or 
to a shell intersection p = r with S, or to the inhnity p = cxd of S. Finally, there is 
a contribution due to the inhnity. 




(cf. Eq.(46) of II). 

In this paper, we shall collect all contributions at each particular boundary and 
transform it in several steps to a general simple form. We shall denote the contri¬ 
bution from the center by ©o, from any shell intersection point r by ©r and from 
the inhnity by ©oo- Thus, we have 

©r = ©0 + X/ ’ 

r 

where 


©0 = -ifdU + gdV + hidd)p=o (6) 

and 

©oo = \im.{fdU + gdV + hido'') — MooEoodt . (7) 

p=oo 

The shortest way to calculate ©r found as yet uses the DNEF gauge; this is, however, 
a singular gauge (cf. Sec. 3 of I), and some justihcation is in order. 

In Sec 3.2.2 of II, gauge transformation of the functions F, G and Hi in the for¬ 
mulae (|l|)-(|D have been calculated starting from the requirement that the Liouville 
form remains invariant. The invariance is in fact pointwise, that is, the integrand of 
the Liouville form in any volume part is invariant at any point p. Calculating the 
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form in a singular gauge such as DNEF coordinates then makes sense. At the points 
where the gauge is regular the integrand has the same value as for a gauge. The 
DNEF gauge is singular where the embedding intersects a horizon; the value of the 
integrand at such a point can be dehned as the limit from the left or from the right 
because the integrand is continuous in a regular gauge. 

This argument shows that the difference of the boundary terms obtained from 
the integration over any volume part is gauge invariant and can be calculated in 
the DNEF gauge. The only point at which caution is necessary is the expression of 
the shell variables r and p in terms of embedding variables and Dirac observables 
because the embeddings are not continuous at the shell. Let us, therefore, generalize 
formulae (30) and (31) of II, 

Pout = -i?(r)Ar(i?,f/)t/'(r) , (8) 

Pin = /2(r)Ar(i?y)E'(r) . (9) 


for the momentum and calculate the corresponding formulae for r. Again, the idea 
is that the value of p and r is gauge invariant and we just have to express it in the 
singular gauge. 

Let us start with pout- Let us label the volume parts of S adjacent to the shell 
by the index K = I, r, I meaning left and r meaning right from the shell. The corre¬ 
sponding pieces of Schwarzschild spacetimes are denoted hj M-Ki the Schwarzschild 
masses by Mk and the maximal extensions oi M.k- We denote the singular 

gauge within Mk by Uk and Vk- The coordinates U and V that occur in Eqs. (H) 
and (|) represent a gauge (see Sec. 2.2 of II). Hence, we have 


Pout 


-R(r) 


mj 



U'(r) 



= -^(r) 


dR dUr\ 


\dUr dp 


p=r 


[dUi dp 


and the generalized equation reads 


Similarly, 


Pout = -Rir)Ar{R,uU') . 


Pin = R{r)Ar{R,vV') . 


( 10 ) 

( 11 ) 


As for r, we have to express it in terms of the coordinates of the shell with respect 
to the singular gauge. Thus, we have immediately 

UK{r) = 

VK(r) = 


4 


Uk , 
Vk , 


( 12 ) 

(13) 










where uk and vk, K = r,l are the new variables to describe the position of the shell. 
As the functions Uk{p) and Vk{p) are monotonous, the relation between r and any 
of Uk and vk is well dehned for each K = I, r. Of course, the four parameters ui, vi, 
Ur and Vr are redundant for the description of the shell position that is determined 
by just one parameter r, but the Liouville form which results at the end of the 
calculation contains only one combination of the four parameters. 

If the shell is marginally bound, that is, if it lies at a Schwarzschild horizon of 
both spacetimes left and right, then either uk or vk diverges. However, since the 
Liouville form is gauge invariant, we obtain the “proper” value for it at these points 
of phase space as follows. First, we calculate the Liouville form in a singular gauge 
for all cases but for the marginally bound ones. Second, we transform it into some 
regular coordinates on the phase space. Third, we take the limits to the marginally 
bound cases. This will be done in Sec. O. 


2.1 Contribution from a single shell 

Let us calculate ©r in the case that p, p = r, is an intersection point of S with a 
single shell. We assume that the shell does not he at the Schwarzschild horizons 
of M-k for K = I, r. Let the embedding dehning S be described by the pair of 
functions {Uk{p),Vk{p)) in Aix for each K. Then all contributions to ©r are (cf. 
Eqs. (51H55) of II): 


©r = 


-Ar 
-Ar 
-Ar 
-Ar 
+ Ar 


fdU + gdV + hido^ + pdr 


pdr = 


RRu + {dU + U'dr 
I In (- ^) RRy + g] (dV + V'dv 




.2 , , 

iln RR^i + HAdR 


RRyU'dr - RR vV'dr 




+ pdr. 


(14) 


The symbol Ar (that is defined in I, below Eq. (19)) contains the indices I and r 
implicitly. 

Now, we can start simplifying and transforming the right-hand side of Eq. (pH]). 
The hrst step is based on the properties of the function R along the shell; it is 
continuous, and it is dehned, from both sides, by 

R{r) := RK{UK{r),VK{r),R) . (15) 
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Differentiating this eqnation, we obtain the relation 


= RK,UK{^){dUK{r) + D^(r)dr^ 

+ RRyK^j) (dVKiy) + V^{r)dr^ + i?x,i(r)ho* • (16) 


Moreover, taking differentials of Eqs. (^) and ([T^) yields 

d(UK{Y)^ = dUKij) + U'j^{r)dr = dux 

and 

d(VK{r)) = dVK{r) + V^(r)cir = dvx ■ 
Hence, Eq. (|T^ simplihes to 


0r = -Ar 


+Aj. 


- In ( j d(^R^^ + Fdu + Gdv + Hido'' 


4 V 

RR^uU'dr - RRyV'dr 


+ pdr. 


(17) 


(18) 


(19) 


Further transformations depend on whether the shell at p is out- or in-going. Let 
us assume that it is out-going; the procedure for the in-going shells is analogous. 

In the next step, we calculate the value of the jumps that occur in Eq. (plQl). 
The idea is that the embedding S is if the transformation between our singular 
gauge and a gauge were known, one could compute the jumps. However, such a 
transformation is outlined by Lemma 2 of 1. At a point of an out-going shell that does 
not lie at a Schwarzschild horizon (either a crossing of a shell and a horizon or a shell 
lying at a horizon), we can apply the Lemma with the result; The transformation 
of a regular gauge U and V to Ui and V) is 


U = Ui, V = Vi 


( 20 ) 


left from the shell, and to Ur and K*, it is determined by 


U — Ur — Ur Ul 


( 21 ) 


and 


Rr{Ur, Vr{V), Mr) = Rl{ui, V, Mi) , (22) 

where Rk{Uk, Vr, Mk) expresses the radius coordinate R in the spacetime M.r as 
a function of the DNEF coordinates Ur and Vr and the mass K = I, r. 
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Let us further recall that we admit only embeddings or else the variational 
principle in II does not lead to the desired equations of motion (cf. Sec. 2 of II). 
Hence, the derivatives U' and V of the embedding functions U (p) and V(p) must 


be continuous across the shell. By differentiating Eq. with respect to V right 
from the shell we hnd that 


dVriy) _ Rm(uuV,Mi) 


dy Rr,V^(ur,Vr{V),Mr) ' 

Eq. (^31) implies that for all points V (p) right from the shell we have 


(23) 


We learn from Eqs. and (^) that 


Ar 


ly \ 


V'(p) . 

(24) 

om the shell: 

V'ip) . 

(25) 

= fi,,K,(r)V,(r) 

(26) 

t p = r, that 

E'l = 0 . 

(27) 


Moreover, Eq. m gives 


A. 


RRuU^ 


+ P = 0 . 


(28) 


Next, let us consider the discontinuity in the logarithm in Eq. (p!9D which can be 
expressed as 


A. 




= In 


U'rV'l 

V'rU'l, 


(29) 


Since U (p) is continuous at p = r, we hnd from Eqs. (^0[) and (^) that 

U'r{r) = U'i{r) . (30) 


Eqs. (1^) and (|30[) can be used to simplify (^) and express it in terms of the 
derivatives Ry at p = r. 


A. 


In — 


V' 


= In 


'Rr,Vr(r) 

.Ri,Vi(r), 


(31) 
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The derivative Rk,Vk be determined in terms of the radius 

coordinate Rk and the mass Mx if the position of the point p is further specified. 
It can he in different quadrants of the extended Schwarzschild spacetimes M.k- The 
quadrants have been labeled by sign pairs {a, {3) in Sec. 2.1 of I. The values of a and 
(3 also determine which DNEF coordinates are well dehned in the quadrant {a, (3): 
they are f/" and V^. 

Suppose that the intersection p lies in the quadrant {ax, Pk) of the Schwarzschild 
spacetime A4x- We must have = ai because of the argument of matching diver¬ 
gences of Sec. 2.2 in I. Then the dependence of the function Rx on the coordinates 
Ux and Vx is given by Eq. (5) of I; 


R — 2Mxf^ 


axpK exp 


'-axK^ +PxV^^' 


AM. 


K 


(32) 


where k is the Kruskal function (cf. Eq. (6) of I). From the dehnition of k, the 
following equations result 

{k,{x) — = X 

and ^ 

K,'(x) = —-—^ . 

Using these equations to calculate Rxy^^ we obtain 


dRx _Pk(^ 2Mx 

2 1 Rk , 


(33) 


We know that the expression 


^ f 2M\ 

a/3 (l - —) 


is always positive. Since a is the same from both sides, the expression on the right- 
hand side of Eq. (|33D has the same sign for both values of K and 


In 


dR dVf^\ _ \l-2Mr/R\ 

“ \l-2Mi/R\ 


.dVp 


dR ) 


= Ar (In 


^In 1- — ^ = Ar (In |i? — 2M|) 


Thus, using Eqs. (|2^), (^8[) and (0), we obtain 


0r = Ar(ln \R - 2M\) - A^{Fdu + Gdv + HM) ■ (34) 


Formula (|3^) holds at each intersection of an out-going shell with S that is not a 


point of Schwarzschild horizon, R = 2Mx- Indeed, the first term diverges at such 
an intersection. 
















In the next step, we use the solutions for the functions F, G and Hi that have 
been found in Sec. 3.3 of II; in the DNEF gauge, 


Fk = 0 , 

Gk = 0 , 

HK^ = -{axur + PKVi^) ^ , 

K = I, r. Using the formulae (13) cinci (|13|) , WG finci tlis-t 


0^ = -id(F2(r)) Ar(hi \R - 2M\) + -A, 


(^au + fiv')dM 


(36) 

(36) 

(37) 


(38) 


Eq. (1^ and the definition of n (Eq. (6) of I) imply 


'OkUk + Pk'‘Jk 


= R -\- 2Mx la 


R 


Substituting this for Pk^k, we can cast Eq. 


Qr- = An 


au^^dM - d(^R'^^ In 


R-2M 


2Mj^ 

as follows: 

+ RdM + 2M In 


(39) 


R-2M 


2M 


dM 


(40) 


For the last two terms in the bracket, however, the following identity can be easily 
shown to hold: 


= d 


R-2M 
1 


+ RdM + 2M In 


R-2M 


2M 


dM 


^{R^ - AM^) In \R - 2M\ - M^ ln(2M) + ^R^ + ]^MR 


Subtracting the corresponding exact form from ©r leaves us with 


0r = An 


audM 


= ttr Ur dMr — Ctl Ul dM„ 


For an ingoing shell we hnd in an analogous way that 


0n = An 


PvdM 


= Pr Ur dMr — Pi Vl dMi 


(41) 


(42) 


(43) 


These two formulae summarize our treatment of single shells. Observe that the 
singularity in 0r due to shell crossing a horizon at p has disappeared: Eq. (^3|) is 
well dehned everywhere at any shell that does not lie at a horizon. 
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2.2 Contribution from a crossing point of two shells 


Let us next consider the case in which the embedding passes through the crossing 
point p of an in-going with an out-going shell. The form ©r at the crossing point 
p = r is given by modihed Eq. ([T^), where we have to write poutdr-l-pindr instead of 
pdr, and where the operator Ar has a different meaning: it does not denote jumps 
across one but across two shells. This ought to be explained in more details. The 
intersection of the two shells defines four spacetime regions around the crossing point. 
Let us denote them by Af K = l,r,u,d, where the subscripts stand for left, right, 
up, down. The embedding passes from the spacetime region on the left, Jlii, to the 
one on the right. Air, without entering in the upper or lower intermediate regions. 
The operator Ar in thus refers to the jump from Aii to Adr- Let us choose the 
coordinates in each region AIk to be the corresponding DNEF coordinates Uk and 
Vk- The shells are described by Eqs. (0 and (1^ but now K = r, I, u, d. Hence, 
the coordinates uk,vk in Aix satisfy Eqs. (p!^ ) - ([l8|) for all K. 

Since i? as a function on spacetime is continuous even at shell crossings, Eq. (HD 
still holds, now for all four values oi K = l,r,d,u. Hence, Eq. (16) also holds, in 
particular for K = I and r = r, so we can again collect all derivatives of R as it has 
been done in Sec. 0 The result is analogous to Eq. (0): 



ri ( u'\ / 


©r = -Ar 


+Ar 

RR,uU'dr - RRyV'dr 

Pout dr -F Pindr 


(44) 


We can even use Eqs. (pTl) and (pdj) for the momenta with Ar meaning the jump 
across the two shells. Indeed, let us start with a regular gauge so that Lemma 1 of H 
is applicable. On the Lemma, the jump in the H-derivative across an out-going shell 
is continuous along the shell and it is zero across the in-going one. As the jump of 
any quantity across two shells is the sum of jumps across each, we have the required 
property, and we can proceed with the transformation to the singular gauge as in 
the case on a single shell. Hence, we obtain the formula: 


0r 



( U'{r)\ 
[ V’ir)) 


d(^R'^{r)^ + F{r)du + G{r)dv + Hi{r)do^ 


(45) 


Let us foliate a neighborhood of the crossing point by a family of embeddings. 
The foliation is described in A4k by Uk{t,p), Vk{t, p) where the crossing point 
corresponds to r = t, p = r; i.e., Uxii, r) = uk and lA(t, r) = vk - Let the foliation 
intersect the ingoing shell at p = rin(r) and the outgoing shell at p = rout('r) so that 
rin(t) = r = rout(t)- In particular, this means that if r < t then Uk{t, rout('7‘)) = Uk 
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for K = I, d and Vk{t, rin(r)) = vk for K = d,r, while if r > t then Vk{t, rin('7‘)) = 
vk for K = l,u and Uk{t, rout('r)) = uk for K = u,r. 

To calcnlate the jnmp in the logarithm term on the right-hand side of Eq. (^51), 
we again apply Lemma 2 of I. At each point of the fonr “legs’” of the crossing 
(trajectories of the shells ontside p), we can End coordinates and obtain eqnations 
for the jnmps in V' and U' analogons to Eqs. (^6]) and (^Ol). Therefore, taking the 


limit r = t, we obtain equations like (^) and (^) for all four spacetime regions 
around the crossing point: 


Rr,v{r)V'r{r) = Ru,v{r) ^(r) , 

Riy{r)V\{r) = Rd,v{r)V'd{r) , 
R^,u{r)U'u{r) = Ri^uir) U'i{v) , 
Rr^u{r) U'r{r) = RdA^)U'd{r) , 


17V (r) 

= U'uiv) , 

(46) 


= U'dir) , 

(47) 

n(r) 

= V\{v) , 

(48) 

V'rir) 

= V'dir) . 

(49) 


Computing Rk,v and Rk,u from Eq. (^) as for the single shell, we can cast Eqs. 
(^)-(^) in the following convenient form: 


V\{i) _ Pu{R(t) - 2M„) u'dr) 

V'Ur) ~ p^lnir) - 2Mr) ' V„(r) 

V",(r) _ A(fl(r) - 2Md) u',{r) 

iTW “ /J,(fl(r)-2M,) ’ UW) 

V„(r) _ a,(fl(r)-2M,) v'^(r) 

U'i(v) a„(fl(r)-2Af„) ’ V\{v) 

U'd(T) ar(fi(r)-2M,) 

U'Ar) aj(fl(r) - 2Mj) ’ rv(r) 


(50) 

(61) 

(62) 

(63) 


Let us view them as algebraic equations determining the derivatives 17' and V in 
terms of 7?(r), the masses M/, M^, Mu, Md and the various indices a and fd. Then 
it is easy to see that not all of them are independent. Multiplying the equations for 
U' in Eqs. (|^ ) and ( |5^ and using the 17'-equations from (|^ ) and (EP lead to the 
condition 


ai\ 

[r{j) - 2Mi) 


[r{v) - 2Md) 

1 

[R{r) - 2Mu^ 

) ar\ 

[R{r) - 2Mr) 


while an analogous procedure for fr'-equations yields 

ft(fl(r) - 2M,) _ Pd(R(T) - 2Md) 
Pu(R(t) - 2Af„) “ P,[R(t) - 2M,) 


(54) 


(55) 
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Let us now recall the argument of matching divergence of Sec. 2.2 in I which implies 
that not all of the indices a and (3 are independent. In general, two regions on 
opposite sides of an outgoing shell should have the same a, and two regions on 
opposite sides of an ingoing shell should have the same /?. In our case, this means 
that au = o.ri cxi = ad, (3u = I3i and f3r = f3d- Therefore, Eqs. (|5^) and both 
reduce to the Dray-’t Hooft-Redmount condition § and Q (cf. Eq. (7) of I), 

R{v) - 2Mi _ R(r) - 2Md . . 

R{r) - 2Mu ~ R{r) - 2M/ ^ ^ 


determining the Schwarzschild radius R{r) of the crossing point in terms of the 
masses of the four regions. These variables are therefore not all independent. 

Let us now return to the logarithm term in Eq. (|45|) and use Eqs. (^ ) -(|53| ) to 
calculate its discontinuity from M.i to Air- The discontinuity 


Ar 


In 



In 


/ f/V(r)n(r) \ 

\V'r{v)U'i{v)) 


(57) 


can be expressed in various ways by using different pairs of equations from the set 
(^)- (|5^) . Of course, because of the Dray-’t Hooft-Redmount condition (|5^) all 
these alternative expressions are equivalent. 

The simplest and most symmetric of these is 


A. 



= In 

R(r) - 2Md 

[ V V'J\ 


R{r) - 2Mu 


(58) 


which brings Eq. (^) to the form 

Or = -^d(R2(r)) In \R{v) - 2Md\ + In |R(r) - 2M„| 

-h ^(arUr + (3rVr) dMr - ^ Ad) dMi . (59) 

Notice that we have again chosen the particular solutions (|35D - (|37D for the functions 
F, G and Hi. 

The last two terms in Eq. (|59| ) can be written with the help of Eq. (|39D. Then, 
applying the identity (41), we obtain our hnal expression for the contribution to the 
Liouville form from a crossing point, involving the discontinuities Ar from the left 
to the right region and At from the lower to the upper region: 


0r = 


-At 


(^ — au + (3v^ dM + - Ar {au + (3v^ dM 


1 1 
= 2 ( “ + fduVn) dMu “ 2 ( ~ ^ 

+ ]-{arUr + IdrVr) dMr - ^ -f- dMi . 


( 60 ) 


Again, this formula holds only if no shell at the crossing is marginally bound. 
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2.3 The general formula for all contributions 


Before we conclude this section let us verify that the above formula for a crossing 
point agrees with the previous formulae and (]4^) for ingoing and outgoing shells. 
Indeed, one can consider an embedding starting from the left and then passing either 
above or below or through a crossing point. If it passes above, it crosses hrst the 
ingoing shell and then the outgoing shell. The total contribution from both shells is 
the sum of the contributions (|4^) and (^31): 


^above [^u Pi Vl ci]\^l T CXr Uj- Uy^ dlVly . 


If it passes below, it crosses the outgoing shell hrst and the ingoing one second. By 
summing the contributions (|4^) and we now get 


^below d]\/[(j^ Oti Xll dAf; T P^^ dAdj' Pd Vd dA/^ . 

Expressions (^) and (|6^ are equivalent, and their symmetric sum indeed recovers 
the result (|60D obtained by letting the embedding pass through the crossing point. 
We can prove the equivalence as follows. 

The difference ©above — ©below can be written as 


©above ©below XldAlf; T ')(^ddMd T 'X^yd]\dy y,.dAfj. , 


where 


Xk ■— —c^kUk + PkVr ■ 


Eq. (|5^ allows us to express y’s in terms of the four masses Mx, 


Xk = 2d?(r) + 4Af;^ In |i?(r) — 2Mx\ — 4:Mk ln(2Af/f) , 


(63) 


if we use formula (^BD to calculate R{r): 

MiMr - MdMy 

Ml + M, - Mi - M. ■ 


(64) 


In this way ©above — ©below is a well dehned form on the four-dimensional space 
spanned by Mi, Mr, Md and My Observe that, from the point of view of canonical 
transformations, our “old momenta” are Mi, Mr and Md, while the “new momenta” 
are Mi, Mr and My and the form ©above — ©below is the differential of a generating 
function T that depends on the old and new coordinates, if we can show that it is 
an exact form. Let us show that. 

First, we observe that ©above — ©below is invariant with respect to the swaps 


Mr ^—> Ml, Md<—> My . 


( 66 ) 
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Hence, it is sufficient to show that 

dxi ^ dXi ^ dxd ^ dxu 

dMd dMi ’ dMr dMi ’ dM^ dMd ' 


Second, we observe that xk depends on the masses with index different from K 

( 66 ) 

(67) 

( 68 ) 


Hence the above equations can 

be written 

dxi 

dR{r) 

dXd 

dR{r) 

dR{r) 

dMd 

dR{r) 

dMi ’ 

dxi 

dR{r) 

dXr 

dR{r) 

c?i?(r) 

dMr 

c?i?(r) 

dMi ’ 

dXd 

dR{r) 

dXu 

dR{r) 

dR{r) 

dMu 

dR{r) 

dMd ■ 


Eq. (H) yields 


2R{t) 


dxK ^ 
dR{r) R{r) - 2 Mk 


(69) 


and Eq. ( |6lD implies 


and 


dR{r) 

dMi 

dR{r) 


= 2 


= 2 - 


{Mr-Md){Mr-M^) 
[Ml + Mr-Md- M„)2 


dMd {Mi + Mr-Md-M^Y' 

the other formulae can be obtained through the swaps (|6^. We also find that 

XMi-Md){Mi-M^) 


R{t) - 2Mi = -2 


and 


R{r) - 2Md = 2 


Ml + Mr-Md- Mu 

{Md-Mi){Md-Mr) 


Mi + Mr-Md-Mu ’ 

the other two differences resulting by the swaps. It follows that 


dR{r) 

dM, 


i?(r) — 2 Mr 


Ml + Mr — Md — Mu 


and 


R{t) - 2Mu 


dR{r) _ 
dMd ~ Ml + Mr-Md- Mu 


(70) 


(71) 


Substituting Eqs. (|6^ , (|70|) and (|7l[) or their suitable swaps into Eqs. (|66|) - (|68[) , we 
easily show the equality if we use Eq. (|56D once more. 
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Since the form ©above — ©below is exact, its integral along a curve depends only- 
on the end points of the curve. Let us calculate such an integral in the M-space. 
Let the curve start at the origin and consist of four straight segments: (0, 0, 0, 0) —> 
(Mrf, 0,0,0), (Md, 0,0,0) ^ (Md,M„0,0), (Md,M„0,0) ^ (M^, 0) and 

(Mrf, Mu, Mr, 0) —>• {Md, Mu, Mr, Ml). An elementary integration yields 


T = -2MiMr + 2MdMu 

_ 2^2 , \Mi - Md\\Mi - Mu\ _ 2 , \Mr - Md\\Mr - Mu\ 

' Mi{Mi + Mr - Md - Mu) " Mr{Mi + Mr - Md - Mu) 

, 2m 2 In \M^-Mi\\Md-Mr\ |M. - M,||M. - M.| 

Md{Mi + Mr-Md- Mu) “ Mu{Mi + M, - M, - ' 


(72) 


Strictly speaking, this formula holds only for that part of the M-space where the 
masses Md and Mu are larger that Mi and Mr because the integration curve would 
otherwise cross the corresponding singularities of the integrand. But one guesses 
that this formula holds in all generic subcases as it stands. We can check this as 
follows. 

The function JF is to be a generating function of the canonical transformation that 
brings the Liouville form ©below into ©above- It depends on the old {Mi, Md, Mr) and 
the new {Mi, Mu, Mr) momenta and the usual generating formulae must hold (as 
some new momenta coincide with some old ones, differences of coordinates appear 
instead of the coordinates themselves): 


dJ^ 

dJ^ 

dJ^ 

m 

dT 

MTr 


{^Pd'^d ^d^d) 1 

(73) 


(74) 

{f^ibvi - aiUi) . 

(75) 

{PrhVr — arUr) • 

(76) 


A simple calculation shows that Eqs. (I^)-(P are valid in all subcases with a 7 ^ 0 
and 6 7 ^ 0 for our guessed T as given by Eq. (72). 

The transformation itself can be calculated as follows. Let us introduce the 
following notation 


gi = -aiui 
Pi = -Ml , 
Qi = -pivi 

Pi = -Ml , 


q2 = (3rVr 
P2 = -M, 
Q2 O^r^r 
P 2 = -M, 


O'd'^d (^d'^d , 
P3 = -Md , 
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Q3 ClLu^u T Pu'^u , 

P 3 = -Mu . 










Then ©below = E -Qndpn and ©above = E -QndPn- Eq. for = d gives 


R{r) = -2p3K 


adl3d exp 


El 


and Eq. (p4D yields 

P3 = 


P 1 P 2 - {Pl +P2-P3)P3f^ 

adlSd exp 

f—y 

1 4p3 y 

P3 - P3« 

adl3d exp 

f —1 



Using Eq. (|39D again, we obtain 

Q3 = 2i?(r)-4P3ln 


P(r) + 2 P 3 


-2P, 


where P(r) and P 3 are given by Eqs. ([77|) and (^). Then, 


( 77 ) 


(78) 


(79) 


Pi = Pi , P 2 = P2 

and Eq. (^^ implies 

Qi = qi- 2P(r) + 4pi In 
Q 2 = q 2 - 2 P(r) + 4 p 2 In 


P(r) + 2pi 

-2pi 

P(r) + 2 p 2 

- 2 p 2 


(80) 


(81) 

(82) 


where again P(r) is to be substituted from Eq. dZH). The desired transformation is 
described by Eqs. ([r7| ) -(l8^ ) ; it is a rather involved one. 

Finally, let us make the following observation. Let us view the operator Ar 
as the jump in the positive spacelike direction and the operator At as the jump 
in the positive timelike direction. This orientation of spacetime is hxed by our 
particular definition of DNEF coordinates. Then, we see that along an outgoing shell 
Ar = —At, and along an ingoing shell Ar = At. By substituting these relations into 
expression ( 0 ) we recover our previous expressions (|4^ ) and (^) as special cases. 
Let us therefore summarize all our results of this section by rewriting expressions 
(^), (^) and ( pDD concisely as 


© 


At 


P* dM 


+ Ar 


T dM 


(83) 


where 


and 


D*_ 

tlK • — 


Tk : = 


—axUK + PrVk 
2 

O^RUk + fdKVR 
2 




This holds in general, with p = r corresponding to the point of intersection of the 
shell(s) with the chosen embedding if none of the shells is marginally bounded. 
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2.4 Contributions from center and from infinity 

For the center and infinity, we shall again utilize the freedom in the gauge choice 
and compute in the DNEF coordinates in the spacetime part surrounding the origin, 
which we denote byAto and the part near infinity, which will be denoted by Atoo- 
For simplicity we shall drop the subscripts 0 and oo from U and V in this section. 


2.4.1 At the center 


The contributions to the Liouville from p = 0 are given by Eq. (^, where the 
functions /, g and hi are defined by Eqs. (0)-(0). The solution (|35|) , (|3^) and (|37D 
for the functions A, G, Hi is trivial in the Minkowski part of spacetime. 


F = 0 G = 0 Hi = 0, (84) 

and therefore Eq. (P) reduces to an expression involving only the logarithm terms 
in /, g and hi. 

At the regular center R = 0 (where p = 0) all embeddings have to satisfy the 
boundary condition 

C(0) = -E'(0), (85) 


which guarantees that the embedded hypersurfaces avoid conical singularities. Eq. 
(^) reduces all logarithm terms in Eqs. (0)-(|3) to zero. The functions /, g and hi 
and consequently the Liouville form at p = 0 is therefore trivial: 

00 = 0 . ( 86 ) 


2.4.2 At infinity 

The contributions to the Liouville from p ^ oo are given by Eq. ©• At infinity, we 
restrict the foliation to be parallel to the T = const, where T is the Schwarzschild 
time coordinate in Aloo; more precisely, we assume that 

R{p) ^ p + 0(p-'), T(p) + 0(p-') , (87) 

for the Schwarzschild coordinates T and R along any embedding. We also assume 
that the foliation parameter t satisfies 

f —^ T 

asymptotically so that A/”oo = 1- Finally, we replace E^o by Moo, Moo being the total 
mass of the system. Then we can write 

©oo = lim (fdU + gdV + hido^) — MoodToo . (88) 

p—>oo 
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The functions /, g and hi have always the same functional form, given by Eqs. m- 
and F, G and Hi are given by Eqs. ( 0 )-( 0 ). where is replaced by M^o- 
Let us start from the only non-trivial contribution, Eq. ( 0 . and use the fact 
that V + U = 2T. As p approaches inhnity the time coordinate T approaches the 
asymptotic time Too, which means that the total contribution from Hi is just 


HidF' = -ToodMoo . (89) 

1) yields an exact form. The 


This term together with the term —MoodToo in Eq. ( 
remaining contributions to the Liouville form in Eq. (^) are proportional to the 
logarithm terms and can be summarized as follows: 


R 


H' 


©oo = Y I I R.udU + Rydv + R,ido 


(90) 


Recall that in a Schwarzschild region of mass M the Eddington-Finkelstein coor¬ 
dinates are related to the coordinates T and R by 


U = T-R*, V = T + R*, 

W ^R + 2Mlii " 


2M 


(91) 

(92) 


As p —> cx), it follows from Eq. (|9^ that R* has the asymptotic expansion 

fi-(p)^p + 2Aflti(^)+0(p-‘). (93) 

Eq. (|9l|) determines how the Eddington-Finkelstein coordinates approach inhnity, 
C/(p) ^ -p - 2Mlti (^) + + 0(p-'), (94) 

r(p) — ^ +P + 2A/ln (^^) + Too + 0(p ^), 
and hence how their derivatives U' and V approach it; 

u'(p) ^ -l + 0(p-'), 

V"(p) ^ 1 + o(p-‘). 


(96) 


(96) 

(97) 


The way in which the differentials dU and dV increase as p —> cx) follows directly 


from Eqs. (H) and (@, 

dU{p) ^ -2dMln + dToo + 2dM + 0{p -^), 

dR(p) ^ +2dMln + dT^ - 2dM + 0{p -^), 


(98) 

(99) 
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while the way in which the logarithm in Eq. (|90D behaves is determined from Eqs. 
(El) and (0): 

In (P) 0{p-^). (100) 


The final terms whose behavior we need to determine are the terms multiplying 
the logarithm in Eq. (]90|). Eqs. (|3^ ) and the fact that we are in the (++) quadrant 
imply that the derivatives of R are related to R and M by 


„ R^^ 1 / 2M\V-U^^ 

R. = — Mi - - 1-- Mi 

’ M ’ 2 V R M 


( 101 ) 


After some simple calculations it follows that 


R,uip) ~p ’ ^y^p) 2 ^ ^^P ’ 


R,i{p) 2Mj - 2Mdn 


( P 


\2M 


+ 


0{p ^In 


( 102 ) 


Putting together Eqs. (ED, (1^) and ( p.02|) it is not difficult to show that 


RfjdU + RydV + R^,(ki\(f>) -* 0(p“‘lii 


(103) 


When Eqs. (|100|) and ( [103|) are used in Eq. (|90|) they yield the final result; the total 
contribution to the Liouville form from infinity vanishes: 


0oo = O. 


(104) 


3 Algebra of Dirac’s observables 

All calculations done and all results obtained as yet in this paper have been entirely 
general in the sense that they hold for a system containing any number of out- and 
in-going null spherical shells. In the present section, we return to our original system 
of two shells. Our final aim is to express the Liouville form in terms of a chosen set 
of Dirac observables for all cases studied in I and so to find the symplectic structure 
of the physical phase space. In the previous section, the Liouville form has been 
reduced to the sum over contributions from shells that intersect a Cauchy surface. 
Each intersection of the surface either with single shells or with crossing points of 
shells contributes according to the general formula (|83|). 
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3.1 The parallel shells: Cases A and B 

The subspacetimes between the shells are denoted hy M-k, where the index K = 
l,m,r (see Fig. 1 in I). The spacetime region M.i is flat and lies entirely in the 
(+) +) quadrant. The masses of the other two regions are Mm and M^. The DNEF 
coordinates of each region A4k are denoted by Uk, Vk- The regions M.m and M.r 
lie in the (+, +) U (—, +) part of their corresponding Schwarzschild extensions. The 
two shells are characterized by the index s = 1,2. The position of the s-shell with 
respect to the chart Uk^ Vk covering the region Af a- is defined by = Vks- By 
applying Eq. (^3]) to case A we And the following total Liouville form: 

Op ^-M^m dflAn T ^r‘2 dAA • 


It is clear from (|105|) that the momentum conjugate to the coordinate Vr 2 is —Mr 
and the momentum conjugate to the difference Vmi — Vm 2 is —Mm- These four 
coordinates describe the physical phase space of the system, which is the subset of 
IR"^ defined by the inequalities 


> 0 , Mr > Mm- 


(106) 


The coordinate Vr 2 can take any values between — cx) and +cx), while Vmi — Vm 2 is 
negative in case A and positive in case A’. 

Different coordinates can be defined by the canonical transformation 


M( 1 ) = Mm , M( 2 ) = Mr -Mm, 

'^( 1 ) = Vmi — Vm2 + Vr2 , '^( 2 ) = Vr2 ■ 

These coordinates are associated more closely to the individual shells: M(i) is the 
mass of the first shell and M( 2 ) is that of the second. Moreover, the coordinate 
n(i) can be considered as the advanced time of the first shell with respect to the 
continuous time coordinate at the past null infinity. That is, we can choose a gauge 
so that Vm 2 = Vr 2 a^d the coordinate V is continuous across the shell. 

In terms of these coordinates the physical phase space is the subset of IR"^ defined 
by 


M(i) > 0 , M(2) > 0 . (107) 

The variables n(2) and n(i) are free to take any values between —cx) and +cx2 but 
"^^( 2 ) > 1 ^( 1 ) in case A while n( 2 ) < n(i) in case A’. In case B and B’, the subspacetimes 
are again denoted by Al/p, K = l,m,r, and the masses of the two non-flat regions 
are Mm and Mr- The left region lies in the (-I-, -|-) quadrant of the flat spacetime 
while the other two in the (-I-, -|-) U (-f, —) part of their Schwarzschild extension. 
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The position of the s-shell is defined by = uks- Using ( |83D we find that the total 
Liouville form is 


(108) 


Op ^mi ^m2 “t“ dTU • 


The description of the phase space is similar to the description above, and analogous 
choices of natural physical coordinates can be made. 

3.2 The crossing shells: Cases C or C’, and subcases 

Cases C and C’ are obtained from each other by interchanging the values 1, 2 of the 
label s. So we can consider only one of these cases. Let us also drop s since the 
positions of the shells are already distinguished by the letters u and v. The four 
spacetime regions are denoted by Aix where K = I, r, u, d. The various subcases are 
pictured in Figs. 4-9 in I. The spacetime region Aii is flat and the masses of the other 
regions are Mr, Mu and M^. The left region lies in the (+, +) quadrant. The other 
regions may he in a union of two or four quadrants of their Schwarzschild extension 
depending on whether each shell lies above, exactly on, or below its corresponding 
horizon. The various possibilities are captured by simple relationships between the 
masses Mr, Mu and M^. For this reason, the indices a and b have been introduced 
in I, dehned by 


a := sgn(Mr - Mu) , b := sgn(Mr - Ma) . 


(109) 


Their meaning is the following: Because of the flatness of the left region the only 
shell boundaries that can he on horizons are the boundaries between Aid and Air 
and between Aiu and Air- If a shell lies above the horizon (corresponding to an 
unbound state) the difference between the Schwarzschild mass on its right and that 
on its left is positive, if it lies below the horizon (bound state) this difference is 
negative, and if it lies exactly on the horizon (marginally bound state) it is zero. 
The indices a and b hence provide all the necessary information and characterize 
the subcases Cab by the nine combinations of their values (+, 0 , —). 

The values of the signs ax and jSx for each case Cab have been found in I: the 
relations 


(y-i — +1 , fdi — +1 , ad — +1 , (du — +1 


hold for all cases; if a 7 ^ 0 , 


ttu = ar = a 
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Let us consider the cases with a 7^ 0: C++, C+o, C+_, C_+, C_o and C _The 

coordinates vi, Vu, Uu and Ur are regular, and we can write the form ©above as follows: 

©above = (“OMu + Vu)dMu + ttUrdM^ . (HO) 

This region of the phase space can, therefore, be considered as a Darboux chart with 
the coordinate pairs 


- auu + Vu, -Mu ; aur, -Mr . (Ill) 

Next, consider the cases with 6 7^ 0: C++, Co+, C_+, C+_, Cq- and C _Here, 

ui, Ud, Vd and Vr are regular, and ©below is 

©below = {Ud - bVd)dMd + bvrdMr . (112) 

Hence, this region of the phase space is covered by the Darboux chart with coordinate 
pairs 


-Ud- bvd, -Md ; bvr, -Mr . (113) 

The two charts overlap in C++UC_+UC^_UC _The (canonical) transformation 


between the two charts in the overlapping region is given by the formulae ([77|), (ffSD, 

(114) 


and (|8^ , where we have to set: 

ga = bVr , q3 = Ud- bvd , P 2 = -Mr , Ps = -Md 


and 


Q2 = aUr , Qs = -aUu + Vu , P2 = -Mr , P3 = -Mu , 
and Pi = Pi = 0. Let us write down the transformation explicitly: 

P2 - Psd 


Q 2 = q 2 + dpsK + 4p2 In 


P2 


Ak 

Qs = -4p3« - V-r(p 2 - P 3 ) In 

K — i 


P2 - P 3 H 


P2 -P3 


P 2 = P2 , 

P 3 = "" 


K — 1 


{P2 - P 3 ) , 


where k is a shorthand for 


K = K 


6 exp 


93 

Jps, 


(115) 


(116) 

(117) 

(118) 
(119) 


( 120 ) 


22 















and the coordinates qi and Qi do not occur in the formulae. The generating function 
T' for the transformation is given by Eq. (72) if we set M; = 0 in it; 


T' 


2 MuMd - 2 M^ In 


\Mr - Md\\Mr - 
Mr{Mr -Md- Mu) 


+ 2 MJ In 


\Mr - Md\ 
Mr-Md- Mu 


+ 2 M^ In 


\Mr-Mu\ 

Mr-Md-Mu' 


The charts cover the whole phase space with the exception of case Coo- 

To hnd coordinates that are regular at Coo, we first use Eq. (| 8 ^) rewritten in 
terms of the three masses, 


^ AMdMu ,,, , {Mr - Md){Mr - Mu) 

Md + Mu-Mr Mr{Md + Mu- Mr) 

and consider the fact that Q 2 diverges at Mu = Mr while q 2 does when Md = Mr- 
It follows that the combination 


Q 2 + 4:Mr In \Mr — Mu 


is always regular. 

Let us introduce a new coordinate q by 

4:MdMu 


Q 2 = q- 


Md + Mu — Mr 


- AMr In 


Mr - Mu 


Md + Mu — Mr 


This should be regular; the motivation for addition of further terms is that 


Vu = q- A{Mr - Mu) In 


Mr-Mu 


( 121 ) 


( 122 ) 


Md + Mu — Mr 

hence, q ^ Vu for Mu —^ Mr- Here, the coordinate is defined by Vu ■= Q 2 + Qs] 
we would have = Vu if we shift Vu and Uu so that Uu = Ur- If we express Q 3 in 
Eq. (|79D in terms of masses. 


Q 3 — 


AMdMu 


+ AMu In 


Mr - Mu 


Md + Mu — Mr 


(123) 


Md + Mu — Mr 

we can see that Eq. (|122|) holds. The choice of Vu instead of q does not, however, 
lead to differentiable components of the Liouville form. 

Let us view Eqs. (p.21|) and 


), together with P 2 = —Mr and P 3 = —Mu as 
transformation to new variables g, Md, Mu and Mr. Then 


0 r = -Q 2 dP 2 - O^dPs 

2Md{Mr - Mu 


= qdMr + 


Md + Mu — M, 


2{M‘^ - M^) 

■d{Mr - Mu) -- —dMd 

^ Md + Mu-Mr 

Mr - Mu 


d 


- 2 (M; - M^) In 


Md + Mu — Mr 
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If we subtract the singular exact form, the rest seems to be regular. But we have also 
to calculate the determinant of the symplectic tensor in order to show the regularity. 
Some simplihcation can be achieved in the coordinates q, x, y and z dehned by 

X = Mr - Mu , y = Mr + Mu , z = Md . 

The corresponding components of the symplectic form = dOr are 


= 1 / 2 , 

^xy 0 : 


= 1/2 , 

^2 _ 

Qrz = 2- - 


yz 


X) 


n 


qz 


Q 


yz 


0, 

^x'^ — xz 

^{z-xy 


and the determinant is 


det 



{x - yf 


^MjMl 

{Md + Mu- Mr)^ ■ 


The C-component Vc of the physical phase space is determined by the boundaries 
(cf. I, Eqs. (9) and (10)) 


Mrf > 0 , > 0 , 


and 


0 < Mr < Md + Mu . 


Hence, is C°° and non-degenerate everywhere in Vq- In particular, case Coo, 
which is dehned by Md = Mu = Mr, is a smooth surface in Vc- 


4 Conclusions and outlook 

In the three papers I, II and the present one, complete sets of Dirac observables 
for the system of two null-matter shells have been found and their Poisson bracket 
have been determined. The result is exceedingly simple and analogous to what 
was found earlier for the single shell |§|, but the calculation itself, in spite of many 
improvements, has still been tedious. We don’t give up hopes that methods exist to 
make the calculation truly simpler. 

Let us turn to the question of quantization. It seems that the similarity of the 
present results to those of |Q suggests that at least some quantization methods can 
be imported from |]^. More specihcally, there does not seem to be any reason why 
each shell of the pair could not bounce at the regular center exactly as the single 
shell did in 0. Some singularities that afflict the classical solutions could so be 
again avoided by the quantum theory. 

On the other hand, the pair of shells is crucially different from the single shell in 
one aspect: the shells can intersect and they have a non-trivial interaction at the 
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intersection. The transformation (|116|) - (|120|) describes the change of the (initial) 
Dirac observables before into the (hnal) ones after the crossing. It can so represent 
the dynamical evolution through the crossing and it ought to be derivable from the 
Hamiltonian of the system. It may be difficult to hnd such a Hamiltonian on account 
of the transformation being rather involved. Perhaps there is another choice of Dirac 
observables that simplihes the transformation. 

This is a technical problem. However, there is another very interesting aspect of 
the crossing that may lead to a great change in results. This is the fact that the 
external shell can become bound after the crossing. For the values of the initial 
observables and such that 


K 


exp 


^3 



the Eqs. ( |118|) and (|119|) yield 
and a black hole forms, while for 


Mr < Mu 


K 




the external shell can reach inhnity. Here, we assume that p^lvz > 1 so Mr > M^ 
and the external shell before the crossing is not bound. This indicates that the pair 
of shells might only sometimes re-expand to inhnity and the other times form a 
“quantum black hole”. Recall that the single shell always re-expands and reach the 
inhnity (cf. [^). The possibility that a black hole forms would, therefore, constitute 
a qualitatively new result. We hope to be able to construct the quantum theory 
soon. 
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